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Abstract. The reaction rates of the key stellar reaction of 14 0(a,p) 17 F have been 
reevaluated. It is thought that the previous l - assignment for the 6.15-MeV state 
is incorrect by a careful reanalysis of the previous experimental data [J. Gomez 
del Campo et al, Phys. Rev. Lett. 86, 43 (2001)]. Most probably, the 6.286- 
MeV state is the key 1~ state and the 6.15-MeV state is a 2~ one, and hence the 
resonance at ^=6.286 MeV (J 7r =l _ ) actually dominates the reaction rates in the 
temperature region of astrophysical interests. The newly calculated reaction rates for 
the 14 0(a,p) 17 F reaction are quite different from the previous ones, for instance, it's 
only about 1/6 of the previous value around 0.4 GK, while it's about 2.4 times larger 
than the previous value around 2 GK. The astrophysical implications have been briefly 
discussed based on the present conclusions. 
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1. Introduction 

Explosive hydrogen and helium burning are thought to be the main source of energy 
generation and a source for the nucleosynthesis of heavier elements in cataclysmic binary 
systems, for example, x-ray bursters, etc. [H El E]. During an x-ray burst (a high 
temperature and high density astrophysical site), Hydrogen and Helium rich material 
from a companion star form an accretion disk around the surface of a neutron star where 
the H and He transferred from the disk begin to pile up. The ap chain is initiated 
through the reaction sequence 14 0(a,p) 17 F(p,7) 18 Ne(a,p) 21 Na [I], and increase the rate 
of energy generation by 2 orders of magnitude [3] . In x-ray burster scenarios, the nucleus 
14 0(ti/2=71 s) forms an important waiting point, and the ignition of the 14 0(a,p) 17 F 
reaction at temperatures ~0.4 GK produces a rapid increase in power and can lead to 
breakout from the hot CNO cycles into the rp-process with the production of medium 
mass proton-rich nuclei [5J |6j [7]. Excepting the 15 0(a,7) 19 Ne reaction, this reaction is 
arguably the most important reaction to be determined for x-ray burster scenarios. 

Wiescher et al. [8] calculated the reaction rates of the 14 0(a,p) 17 F reaction, and 
shown that the resonant reaction rates dominated the total rates above temperature 0.4 
GK. However, Funck et al. [HI [10] found that direct-reaction contributions to the £=1 
partial wave are comparable to or even greater than the resonant contributions at certain 
temperatures. Because the resonant reaction rates of 14 0(a,p) 17 F depend sensitively on 
the excitation energies, spins, and partial and total widths of the relevant resonances in 
18 Ne, Hahn et al. [H] extensively studied the levels in the compound system 18 Ne [TT] 
by several reactions, such as, 16 0( 3 He,n) 18 Ne, 12 C( 12 C, 6 He) 18 Ne as well as 20 Ne(p,t) 18 Ne 
reactions. Based on the firmer experimental results, they concluded that this reaction 
rate, in the important temperature regime ~0.5-l GK, was dominated by reactions on 
a single 1~ resonance at an excitation energy of 6.150 MeV lying 1.036 MeV above the 
14 0+a threshold (Q Q =5.114 MeV). Harss et al. [12J studied the time reverse reaction 
17 F(p,a) 14 in inverse kinematics with 17 F beam at Argonne, and identified three levels 
at 7.16, 7.37, 7.60 MeV and determined their resonance strengths as well. Later, 
Gomez del Campo et al. [13] used the p( 17 F,p) resonant elastic scattering on a thick 
CH 2 target to look for resonances of astrophysical interest in 18 Ne at ORNL. In the 
region investigated, they located four resonances at excitation energies of 4.52, 5.10, 
6.15, and 6.35 MeV in 18 Ne, and ^=1", 2~ were respectively assigned to the last two 
states based on their i?-matrix analysis. Subsequently, Harss et al. [H] extracted the 
resonance strength and the width T a for the 6.15-MeV state based on this 1~ assignment 
together with the excitation function obtained from their previous work [12]. Recently, 
the inelastic component of this key 1~ resonance in the 14 0(a,p) 17 F reaction has been 
studied by a new highly sensitive technique at ISOLDE/CERN [15], and found that this 
inelastic component will enhance the reaction rate, contributing approximating equally 
to the ground-state component of the reaction rate, however not to the relative degree 
suggested in Ref. [16] . 

As a summary, all the previous discussions and calculations [HI [TJJ [161 US] related 
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Figure 1. The upper panel of Fig. 2 as shown in Ref . |13j . The present 7?-matrix fit 
with the same resonant parameters is shown for comparison. 



to the reaction rates of 14 0(a,p) 17 F are based on the 1~ assignment for the 6.15-MeV 
state. In this work, the spin-parities for three relevant states (at 6.15, 6.286, and 6.345 
MeV) will be reassigned by a carefully reanalyzing the experimental data measured at 
ORNL [13], and it is found that most probably the 6. 286- MeV state is the key 1~ state 
and the 6.15-MeV state is a 2~ one. Therefore the resonance at E X =Q. 286 MeV ( J % =1~) 
actually dominates the reaction rates in the temperature region of astrophysical interests 
rather than the 6.15-MeV state does. The reaction rates have been recalculated and 
the astrophysical consequences have been briefly discussed based on the present new 
assignments. 

2. Data Reanalysis 

The previous ORNL experimental data[13], i.e., the integrated cross sections of the 
17 F+p elastic scattering in the angular range of #cm = 162°~178°, have been reanalyzed 
by a R- matrix method [171 [TBI [191 [21] (see example [20]) carefully. A channel radius 
of ro=1.25 fm (R=roX (1+17 1 / 3 )) appropriate for the 17 F+p system [TBI 181 111] has been 
utilized in the present .R-matrix calculation, where the fitting results are insensitive to 
the choice of radius, e.g., the resultant x 2 /N value is changed by less than 15% in the 
region of r =1.2~1.3 fm. The ground state spin-parity configurations of 17 F and the 
proton are 5/2 + and l/2 + , respectively. Thus, there are two channel spins in the elastic 
channel, i.e. s = 2, 3. But for a 1~ state, only one unique channel spin, s—2, is possible 
for an i=\ transfer. 

The fitting curve shown in the upper panel of Fig. 2 in [13] has been exactly 
reproduced by utilizing their resonant parameters as shown in Fig. [TJ i.e., E cm =0.6 MeV, 
3+ T=18 keV; E cm =\.\8 MeV, 2+ T=45 keV (1.18 MeV was mistyped by 1.118 MeV 
in [13]); and the last 'groove' structure can be reproduced by the following parameters: 
£ cm =1.53 MeV, 2", T=5 keV, which is consistent with the results of Ref. [TT](r<20 




Figure 2. i?-matrix fits for the experimental data measured at ORNL [T3] for the 
clastic scattering of 17 F+p. The vertical scale corresponds to the angle integrated cross 
sections in the range #cm=162°^178°. (a) fitting for two resonances, (b) fitting for 
three resonances. All curves are convoluted by an assumed 10-keV energy resolution 
except those two labeled by 'fitl' and 'fit2'. For comparison, the ORNL fit is shown 
as well. See text for details. 



keV). In Fig. [H no energy resolution is convoluted in the fitting curve, and it implies 
that this effect was not considered (or neglected) in the previous analysis [13]. The data 
shown in the lower panel of Fig. 2 in [13j have been carefully analyzed, and it's found 
that the 1~ assignment for the 6.15-MeV state is very unlikely as shown in Fig. 2(a) (see 
curves labeled by 'fitl' and 'fit2') with Gomez del Campo et al.'s parameters listed in 
Table 1. Here no energy resolution is convoluted in the curves labeled by 'fitl' and 'fit2'. 
It shows that the shape of the 1~ resonance (^=1) is of a 'groove' structure instead of a 
'bump' one. In addition, all possible combinations of different spin-parity assignments 
for these two states have been attempted and the most probable fitting curves are shown 
in Fig. 2(a) with parameters listed in Table 1. In order to achieve a better fit, all curves 
are convoluted by an assumed 10-keV energy resolution, of course this will not affect 
the spin-parity assignments for the resonances. Furthermore, we have tried to fit the 
data with three resonances, and the most probable fitting curves are shown in Fig. 2(b) 
with parameters listed in Table 1. It's very obvious this kind of three- resonance fits 
reproduce the experimental data better than those two-resonance ones, especially the 
'groove' structure at E cm =2.32 MeV can be well fitted (see x 2 /N values in Table 1). The 
present i?-matrix analysis shows that two states at E cm =2. 20(E X =6. 12), 2.40(£' x =6.32) 
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Table 1. Resonant parameters used in Figj2j Here, resonance energies (E r ) are in 
units of MeV, and proton partial widths (T p ) in keV. The parameters in 'fitl' and 'fit2' 
were used in the previous work [13J . See text for details. 



Resonance 1 



Resonance 2 



Resonance 3 



Sets x 2 /N 



E rl J*[£,a] T p E r2 J"[£,a] T p E r3 J" [I, a] T p 



fitl a 


31.5 


2.22 


1" 


[1,2] 


50 


2.42 


2" 


[1,3] 


50 










fit2 a 


15.1 


2.22 


1- 


[1,2] 


50 


2.42 


2" 


[1,2] 


50 










fit3 b 


2.0 


2.20 


3- 


[1,3] 


15 


2.39 


3" 


[1,2] 


30 










fit4 b 


1.8 


2.20 


2~ 


[1,2] 


15 


2.41 


T 


[1,2] 


30 










fit5 b 


3.1 


2.20 


3- 


[1,3] 


10 


2.40 


T 


[1,2] 


20 










fit6 b 


2.4 


2.20 


2" 


[1,2] 


20 


2.41 


3" 


[1,3] 


10 










Setl b 


1.2 


2.20 


3- 


[1,3] 


12 


2.40 


2" 


[1,3] 


12 


2.32 


1- 


[1,2] 


15 


Set2 b 


1.1 


2.20 


2- 


[1,2] 


20 


2.40 


2" 


[1,2] 


20 


2.32 


1- 


[1,2] 


15 


Set3 b 


1.2 


2.20 


3- 


[1,3] 


10 


2.41 


3" 


[1,3] 


12 


2.32 


1- 


[1,2] 


10 


Set4 b 


1.1 


2.20 


2- 


[1, 3] 


15 


2.41 


3" 


[1,2] 


10 


2.33 


1- 


[1,2] 


12 



a No energy-resolution convoluted in the fit curves of Figj2] 

b A 10-keV energy-resolution convoluted in the fit curves of Fig|2] 



both possibly have J 7r =2~ or 3~, while the state at E cm =2.32(E x =Q.2A) most probably 
has J 7r =l _ . These three states should correspond to the E x =6.150, 6.345, and 6.286 
MeV states observed before [IT] within a ~30 keV uncertainty 

3. Discussion 

In order to constrain the spin-parity assignments for these states, let's examine the 
well-known mirror nucleus 18 0, which has only three known levels J 7r =l~,(2~), and 3~ 
around this energy region [II] . According to the previous results [15] [16], the 6.15 MeV 
state can decay to the first excited state (l/2 + ) in 17 F with an appreciable width. In this 
case, two channel spins are s=0, 1, respectively. Such a decay is possible only for £=3 
but, the barrier penetrability is extremely too small to produce such a large decay width. 
Thus, this 6.15 MeV state most probably is a 2~ state since the 1~ assignment has been 
excluded in the present i?-matrix analysis. The results from the 12 C( 12 C, 6 He) 18 Ne and 
20 Ne(p,t) 18 Ne reactions suggest that 6.286-MeV state is of natural parity and 6.345-MeV 
state of unnatural parity. Therefore, we propose that these two states most probably 
have 1~ and 2~, respectively Actually Funck et al. [9] predicted a 1~ state at 6.294 
MeV. The ORNL experimental data can be reproduced very well with these assignments 
(see fitting curve labeled by 'Set2' in Fig. 2(b)). 

Therefore only one natural-parity state, i.e., ^.=6.286 (J""=l~, £ a =l), is needed 
to calculate the reaction rates of 14 0(a,p) 17 F in the temperature region interesting for 
x-ray burster scenarios [HI [HI [15], Its r Q partial width is about 13.4 eV according 
to the relationship of r a oc C 2 S a x Pe(E r ) [22], while it was 0.34 eV in the previous 
work [H]; its resonant strength W7( a|P ) is about 40 eV rather than the previous 2.4 eV. 
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Figure 3. The ratios between the present resonant reaction rates and those previous 
ones at certain temperature range. See text for details. 

The resonant reaction-rate ratios between the present results and the previous 
ones [H] are plotted in Fig. |3j Only one resonance (E x =6.286 MeV) is included in 
the present work while two resonances (E X =Q.150, 6.286 MeV) involving in the previous 
work. It can be seen that the present reaction rate is quite different from the previous 
one, for instance, it's only about 1/6 of the previous value around 0.4 GK but, it's about 
2.4 times larger than the previous value around 2 GK. According to the present analysis, 
we think the 1~ assignment for the 6.286-MeV state is reasonable, and hence the spin- 
parity assignments for the 6.150, 6.345 MeV states are rather unimportant in calculating 
the reaction rates (i.e., whether they are J 7r =2~,3~, or vice versa). The present rates 
confirms that the 14 0(a,p) 17 F reaction is rather unlikely to be dominant component 
in the hot CNO cycles in novae environments (instead 15 0(a,7) 19 Ne reaction is, see 
discussions in [H]). Due to the present rate enhancements above 0.9 GK, this reaction 
can, however, contribute strongly to the breakout from the hot CNO cycle under the 
more extreme conditions in x-ray bursters. The present conclusion could probably affect 
the onset temperature where the a-capture dominates /3-decay and a breakout from the 
hot CNO cycle via 14 0(a,p) 17 F reaction begins to take place jH]. However, it should 
be noted that the statistics of ORNL data is not very good, and hence it prevents us 
from putting the present assignments on a very firm ground. Therefore more precise 
experiment, e.g., (p,t) transfer reaction, is strongly suggested to confirm these results 
in the future. 
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